This paper presents detailed analysis and modelisation of 
the star HD 17051 (alias t Hor), which appears as a specially 
interesting case among exoplanet host stars. As most of these 
stars, t Hor presents a metallicity excess which has been mea- 
sured by various observers who give different results, ranging 
from [Fe/H] = 0.11 to 0.26, associated with different atmo- 
spheric parameters. Meanwhile the luminosity of the star may 
be determined owing to Hipparcos parallax. Although in the 
southern hemisphere, this star belongs to the Hyades stream 
and its external parameters show that it could even be one of the 
Hyades stars ejected during cluster formation. The aim of this 
work was to gather and analyse our present knowledge on this 
star and to prepare seismic tests for future observations with 

■ the HARPS spectrometer (planned for November 2006). We 
(^) have computed evolutionary tracks with various metallicities, 

in the two frameworks of primordial overmetallicity and accre- 
\ tion. We have concentrated on models inside the error boxes 

■ given by the various observers in the log g - log T e gf diagram. 
/"■n ' We then computed the adiabatic oscillation frequencies of these 

■ models to prepare future observations. The detailed analysis 
of t Hor presented in this paper already allowed to constrain 

,— ' its external parameters, mass and age. Some values given in 
' the literature could be rejected as inconsistent with the over- 

■ all analysis. We found that a model computed with the Hyades 
t— I ' parameters (age, metallicity) was clearly acceptable, but other 

• ones were possible too. We are confident that observations with 
HARPS will allow for a clear conclusion about this star and that 
it will bring important new light on the physics of exoplanet 
\Q ' host stars. 
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Abstract. 

1. Introduction 

Studying the internal structure of exoplanet host stars is partic- 
ularly important for the understanding of planetary formation. 
In this framework, asteroseismic studies represent an excellent 
tool to help determining the structural differences between stars 
with and without detected planets. Among these differences, 
the observed overmetallicity of exoplanet host stars compared 
to other stars (Santos et al. |2"U0H and l2005l Gonzalez |2003 
Fischer and Valenti 2005 1 needs to be understood. Two extreme 
scenarios are still possible to account for this observed over- 
metallicity : the primordial origin, which assumes that the stars 
formed out of an overmetallic nebula, and the accretion origin 
for which the observed metallicity is due to accretion of hydro- 
gen poor material onto the star during planetary formation (see 
Bazot and Vauclair 2004 1. 

These authors pointed out that the evolution of stars with 
masses around 1 . 1 M is very sensitive to their internal metal- 
licity, due to the possible formation of a convective core. In this 
particular region of the HR diagram, main- sequence stars with 
solar internal metallicities have no convective cores while over- 
metallic stars develop convection in their central regions: evo- 
lutionary tracks computed with different internal metallicities 
may cross the same point in the HR diagram while they cor- 
respond to models of quite different masses and different past 
histories. This behavior was used to try to derive whether the 
exoplanet host star /i Arae is overmetallic from its surface down 
to its center (overmetallic scenario) or only in its outer layer 
(accretion scenario). The star was observed with the HARPS 
spectrometer in June 2004. Up to 43 p-modes could be iden- 
tified (Bouchy et al. 12005 1 and a detailed modelisation could 
be achieved (Bazot et al. I2005> . A possible test to determine 
the internal metallicity of the star, using the frequency small 
separations, was discussed. Unfortunately the modes with fre- 
quencies around and above 2.5 mHz could not be identified 
with enough precision to reach a definitive answer. More so- 
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Table 1. Effective temperatures, gravities and metal abun- 
dances observed for t Hor. [Fe/H] ratios are given in dex. The 
references for the given values are given in column 4. 



T effW 


logg 


[Fe/H] 


Reference 


6136±34 


4.47+0.05 


0.19±0.03 


Gonzalez et al. 2001 


6252±53 


4.61±0.16 


0.26±0.06 


Santos et al. 2004 


6097±44 


4.34±0.06 


0.11±0.03 


Fischer & Valenti 2005 



phisticated signal processing techniques should be developped 
to try to go further in this direction for fi Arae. 

As will be discussed in the present paper, the situation for 
the star t Hor (HD 17051, HR 810) is quite different from that 
of ju Arae. With a slightly larger effective temperature, suggest- 
ing a larger mass (between 1.14 and 1.22 M G , as will be seen 
below), this star may be able to develop a convective core for a 
solar internal metallicity as well as for larger ones. The evolu- 
tionary tracks corresponding to the overmetallic and accretion 
scenarios both show the characteristic behavior induced by the 
presence of a convective core during evolution. 

However, among planetary host stars, i Hor is a special 
case for several reasons. Three different groups have given dif- 
ferent stellar parameters for this star, summarized in table 1. 
Meanwhile, Santos et al. 2004 suggest a mass of 1.32 M for 
this star while Fischer and Valenti 2005 give 1.17 M . 

This star belongs to the "Hyades stream", which means that 
it has the same kinematical characteristics as the Hyades clus- 
ter in the Galaxy, although it is located several tens of parsecs 
away from the cluster itself (Chereul et al l 19991 Grenon^OOO 
Chereul and Grenon 2000 Kalas and Delorn 2006 1. There are 
two possibilities to explain this behavior : either the star has 
been captured inside the stream by tidal effects in the Galaxy, 
or it can have been formed together with the cluster, and evapo- 
rated during cluster formation. In this case, the star should have 
the same age and internal metallicity as the Hyades stars. As 
will be seen below, asteroseismology can solve this question, 
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which can give important clues for the general understanding 
of exoplanets-host stars. 

For the present paper we computed evolutionary tracks and 
models lying inside the error boxes as given by the observers 
for this star, in the log g - log T e ff and log L/L Q - log T e fj 
planes, using the TGEC code (Toulouse-Geneva Evolutionary 
Code) as described in previous papers (e.g. Bazot and Vauclair 
2004 and Bazot et al. l2005> . We tested the differences obtained 
in the internal structure of the models for the various possi- 
ble observed metallicities and the overmetallic versus accretion 
scenarios. We computed the oscillation frequencies for sev- 
eral characteristic models with the aim to discuss asteroseismic 
tests and suggest specific observations. 

The present computations already lead to important results 
concerning the metallicity, mass and age of this star. We find 
that a metallicity larger than [Fe/H] = 0.20 dex is quite unprob- 
able and that the stellar mass cannot exceed 1.22 M , contrary 
to Santos et al. 2004 estimate. Furthermore, models with the 
metallicity and the age of the Hyades are quite realistic. We 
show how asteroseismic observations of this star could solve 
this question: if we can observe this star and identify the oscil- 
lation mode frequencies, we should be able to derive whether 
this star has been formed together with the Hyades or not. 
While a negative answer would be interesting in itself, a posi- 
tive answer would be taken as a proof that the overmetallicity 
has a primordial origin. 

The computations of evolutionary tracks and various mod- 
elisations of this star are described in section 2. Some conclu- 
sions are already given at this point. Section 3 is devoted to the 
model computations, the asteroseismic tests and predictions. A 
summary and discussion of the results are given in section 5. 

2. Evolutionary tracks and models 

2.1. Observational boxes and computations 

Three different groups of observers have determined the metal- 
licity and external parameters (T e gp, log g) of l Hor : Gonzalez 
et al. l2"U0Tl Santos et al. I2U041 and Fischer & Valenti l2005l (see 
tableQ). 

Gonzalez et al. 2001 have observed t Hor with the CTIO 
1.5 m with the fiber-fed echelle spectrograph. To compute the 
stellar parameters (T e fj, log g and [Fe/H]), they use the line 
analysis code MOOG (Sneden fT973l . the Kurucz (1T9931 LTE 
plane parallel atmospheres, and Fe I and Fe II equivalent width 
(EW) measurements. The gf-values of the iron lines were cal- 
culated from an inverted solar analysis using the Kurucz et al. 
( 1984 1 Solar Flux Atlas and their spectrum of Vesta. 

Santos et al. |2004 have obtained spectra with the FEROS 
spectrograph (1.5-m and 2.2-m ESO/MPI telescope, La Silla, 
Chile). They have derived stellar parameters using a standard 
LTE analysis with the code MOOG (Sneden 1 973k a grid of 
Kurucz dT9^3l ATLAS atmospheres and 39 Fe I and 12 Fe II 
lines. The log gf values for the iron line were computed from an 
inverted solar analysis using the Kurucz et al. ( 1984 1 Solar Flux 
Atlas and a Kurucz grid models for the sun (Kurucz 1993 1. 

The spectra used by Fischer & Valenti 2005 were ob- 
tained with the Anglo-Australian Telescope. The computations 



of stellar parameters were done assuming LTE. They created 
a synthetic spectrum with a radiative transfer code using the 
Kurucz stellar atmosphere models (Kurucz Tl993t and with an 
atomic line data basis (Vienna Atomic Line Data-base [ VALD] ; 
Kupka et al 1999 Ryabchikova et al. 1999}. This code uses a 
fitting algorithm to obtain T e fj, log g, v sin i and abundances. 
They have made small adjustements to obtain astrophysical 
atomic line parameters (log gf values) by fitting to the Kurucz 
et al. ( 1984 1 Solar Flux Atlas and by analysing several Vesta 
observations. 

This star has also been observed with the Hipparcos satel- 
lite from which the parallax was derived: n = 58.00 + 0.55 mas. 
The visual magnitude of i Hor is given as V=5.40 (SIMBAD 
Astronomical data base). The overall interval of effective tem- 
peratures obtained from the three observing groups is: T e g = 
6179 ± 126. Using this as an uncertainty on temperatures and 
with the tables of Flower (Flower 1996 1, we obtained for the 
bolometric correction: BC= -0.023 ± 0.01. With a solar ab- 
solute magnitude is M/, i,e = 4.746 (Lejeune et al. 1998), we 
deduce a luminosity of log L/L = 0.219 ± 0.024 

We computed series of "overmetallic" and "accretion" 
models which could account for the observed parameters of t 
Hor. We used the Toulouse-Geneva stellar evolution code, with 
the OPAL equation of state and opacities (Rogers & Nayfonov 
2002, Iglesias & Rogers 1996) and the NACRE nuclear re- 
action rates (Angulo et al. 1999). In all our models micro- 
scopic diffusion was included using the Paquette prescription 
(Paquette et al. 1986, Richard et al. 2004). The treatment of 
convection was done in the framework of the mixing length 
theory and the mixing length parameter was adjusted as in the 
Sun (a = 1.8). The effect of changing this value is discussed 
in section 2.4. For the overmetallic models, we considered two 
different cases for the helium value. In the first case, we as- 
sumed that helium was enriched as well as metals according 
to the law given by Isotov & Thuan 120041 dY/dZ =2.8 + 0.5. 
In the second case, we assumed that the primordial cloud was 
only metal enriched, with a solar helium value. The accretion 
models were computed with the same assumptions as in Bazot 
and Vauclair 2004 (instantaneous fall of matter at the beginning 
of the main sequence and instantaneous mixing inside the con- 
vection zone). Neither extra-mixing nor overshoot were taken 
in account in the present paper. 

2.2. Evolutionary tracks for overmetallic models 

We first present the evolutionary tracks computed for over- 
metallic models, that is models with initial overmetallicities, 
using the three different values given by the observers : [Fe/H] 
= 0.11, 0.19 and 0.26. In these models, the mixing length pa- 
rameter is a = 1.8, adjusted as in the Sun. 

Figure 1 displays the results obtained in the log L/L - 
log T e fj frame (on the left) and in the log g - log T e fj frame (on 
the right). In each graph, the three error boxes are drawn, but 
the one which corresponds to the chosen metallicity (as given 
by the same authors) is presented in thicker lines. Tracks com- 
puted for enhanced helium are drawn in solid lines while those 
computed with solar helium are represented by dotted lines. 
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Some differences can be observed between the two presen- 
tations (log L/L - log T e ff and log g - log T e ff planes). This 
is due to the fact that the error boxes do not correspond to 
the same observable parameters. The luminosity is computed 
from the Hipparcos parallax and from the visual magnitude: it 
does not depend on the other observed parameters. On the other 
hand, the log g values are directly related to the other observed 
parameters as determined by each observing group. In the fol- 
lowing, we choose to work in the log g - log T e ff plane, which 
is more consistent than the log L/L - log T e ff plane from the 
observing viewpoint. 

We note from Figure 1 that none of the evolutionary tracks 
cross the Santos et al 2004 error box. We will come back on 
this result below (section 2.4). The masses of the models which 
may be found in error boxes are in the range 1 . 1 2M to 1 . 1 8M 
for enhanced helium, and in the range 1.1 4M to 1.22M for 
solar helium. 

2.3. Evolutionary tracks for accretion models 

In Figure 2 (left), the evolutionary tracks of models with accre- 
tion are displayed for the same three metallicities. The masses 
of the models consistent with the observed parameters lie be- 
tween 1.1 1M and 1.1 6M : they are smaller than those of the 
overmetallic models, while the ages are similar. 

For such stellar masses, a convective core appears during 
stellar evolution in the accretion models as well as in the over- 
metallic models, which explains that both kinds of tracks have 
the same characteristic behavior. We will see however in sec- 
tion 3 that the accretion models which could account for t Hor 
are too young for the convective core to be developped in the 
accretion models, while it already exists in the overmetallic 
ones 

We note that for the accretion models as well as for the 
overmetallic ones, no solution can be found in the Santos et 
al. 2004 error box, which correponds to a high metallicity 
[Fe/H] = 0.26: the metallicity of i Hor is certainely lower. 

2.4. The effect of changing the mixing length 
parameter 

The evolutionary tracks displayed on Figures 1 and 2 (left) have 
been computed with a solar mixing length parameter (a — 1 .8). 
We studied the effect of increasing the value of this parameter. 
Figure 2 (right) displays overmetallic models computed with 
a = 2.2. For the same position in the log g - log T e gp diagram, 
the models are more evolved and less massive than for a smaller 
mixing length parameter. In this case, we can find very young 
models in the Santos et al. 2004 error box. 

Such a large mixing length parameter does not seem realis- 
tic for this star. Israelian et al. |2t)04 have derived its lithium 
abundance: they give loge(Li)=2.63, a somewhat high value 
which is of the order of the value determined for the Hyades in 
this spectral range. It is also consistent with the computations 
of lithium depletion obtained with a mixing length parameter 
adjsuted on the Sun. A larger lithium depletion would be ex- 
pected with an increased depth of the outer convection zone. 



However, as these results still depend on the details of hydro- 
dynamics, with unknown parameters, we did not eliminate the 
models computed with a large a and kept one of them for as- 
teroseismic studies (model OM3, Table 2). 

3. Models 

3.1. Choice of models 

We chose six different models which all could account for the 
observable parameters of l Hor, to prepare their seismic analy- 
sis (Tables 2 and 3). Overmetallic models OM1 and OM2 and 
accretion models AC1 and AC2 lie in the Fischer and Valentini 
2005 and Gonzalez et al . 1200 II error boxes. Models OM2 and 
AC2 have a luminosity slightly too high compared to that deter- 
mined from the Hipparcos parallax (log L/L = 0.219 + 0.024). 
It was not possible to find models satisfying both constraints of 
log g and log L/L for this metallicity. Model OM3 has been 
computed with a larger mixing length parameter a — 2.2 and 
lies in the Santos et al 2004;rror box. Finally model OM4 is 
a special overmetallic model which has been computed to fit 
as precisely as possible the parameters of the Hyades: we have 
chosen an age of 627 Myr (Perryman et al. 1998 land a metal- 
licity [Fe/H]=0.14 (Cayrel de Strobel et al. [T997l . 

3.2. Seismic studies 

Adiabatic oscillation frequencies were computed for the mod- 
els chosen in section 2, using the PULSE code (Brassard et al. 
1992 The frequencies were computed for angular degrees 1=0 
to 1=3 and radial orders ranging typically from 4 to 100. 

We first computed the large separations, Av = v n+ ij - v n -\,i 
(Table 2), and drew the echelle diagrams for the six selected 
models (Figure 3). We then computed the small separations, 
5v - v„j - v„-ij + 2, which are presented in Figure 4. All 
these computations are similar to those presented in Bazot and 
Vauclair I2"0"0l1 and Bazot et al l2005l for n Arae. 

The star yu Arae was observed during height nights with 
the HARPS spectrometer in La Silla in June 2004. The pre- 
cision of the observations was such that the large separation 
could be determined with an accuracy better than 1 yuHz and 
the observed echelle diagram could be rpecisely compared with 
thise obtained from modelisation. We are confident that for l 
Hor, similar observations will allow an unambiguous selecion 
among the models presented here. 

4. Summary and discussion 

We presented in this paper a detailed modelisation of the 
exoplanet-host star ( Hor, in view of future seismic analysis. 
As most of these stars, i Hor has a larger metallicity than stars 
without planets. A special interest for this star is due to the fact 
that it belongs to the Hyades stream and could possibly have 
been formed together with the cluster stars and ejected during 
their formation process. 

We computed models with the two extreme hypothesis of 
primordial overmetallicity on the one hand (overmetallic mod- 
els), accretion induced overmetallicity on the second hand (ac- 
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Fig. 1. These graphes display the error boxes for the position of the star l Hor in the log L/L Q - log T e ff diagram (left) and in the 
log g - log T e ff diagram (right), as given by the three groups who have observed this star (Table 1). The luminosity is obtained 
from the Hipparcos parallax (see text). Evolutionary tracks are computed for completely overmetallic cases. Those computed 
with enhanced helium are drawn in solid lines, while those computed with solar helium appear in dashed lines. The three error 
boxes are: Santos et al. 2004 (diamonds), Gonzalez et al. 2001 (triangles), Fischer & Valentini 2005 (asterix). The author's boxes 
which correspond to the chosen metallicity, in each graph, are enhanced in boldface lines. 



cretion models). We gathered the various determinations of the 
observable parameters for this star, as given by several observ- 
ing groups : Gonzalez et al. 2001 Santos et al. 2004 and Fischer 
and Valenti 2005 The luminosity of this star was obtained us- 
ing Hipparcos parallax. 

In this first step of modelisation, we were already able to 
constrain the values of the external parameters, compared to 
those given by the observers, for a simple reason of consistency 
of the whole set of parameters. A high metallicity of [Fe/H] = 



0.26, as given by Santos et al. 2004 seems excluded. We expect 
a value between 0. 1 1 , as given by Fischer and Valenti 2003and 
0.19, as given by Gonzalez et al. 120011 We can also exclude 
masses larger than 1.22 M : we find that the mass of this star 
should lie between 1.14 and 1.22 M . 

We have computed the oscillation frequencies, the large 
separations, the small separations and we have drawn the 
echelle diagrams for six possible models of l Hor. Two of these 
models (OM1 and OM2) are overmetallic, with metallicities 
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Fig. 2. Evolutionary tracks in the log g - log T e ff diagram for the three metallicities, as in Figure 1 . Here are presented evolutionary 
tracks obtained with the accretion hypothesis (left) and evolutionary tracks obtained with an increased mixing length parameter: 
a = 2.2 instead of a = 1.8. The error boxes are the same ones as in Figure 1. 

Table 2. Mass, age, gravity, effective temperature, luminosity, surface metallicity, acoustic depth and large separations for stellar 
models satisfying to the observational constraints for i Hor. 



Model 


M* 


(Mo) 


Age (Gyr) 


logg 


logT eff (tf) 


log L/L 


[Fe/H] 


tac (s) 


Ay QMz) 


OM1 


1 


140 


2.231 


4.38 


3.786 


0.207 


0.11 


4167 


120 


OM2 


1 


180 


0.522 


4.43 


3.787 


0.185 


0.19 


3876 


129 


AC1 


1 


120 


2.979 


4.36 


3.786 


0.225 


0.11 


4310 


116 


AC2 


1 


150 


0.522 


4.43 


3.788 


0.170 


0.19 


3817 


131 


OM3 {a = 2.2) 


1 


190 


0.508 


4.46 


3.797 


0.191 


0.26 


3650 


137 


OM4 (Hyades) 


1 


180 


0.627 


4.42 


3.790 


0.202 


0.14 


3906 


128 
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Table 3. Radii, depth and mass of the convective core, radius at the bottom of the outer convective zone, mass of the outer 
convective zone, initial and present surface chemical composition for stellar models satisfying to the observational constraints 
for t Hor. 



Model 


R* (cm) 


WR* 


M cc /M* 


r ec/R* 


M ec /M 4 


Y 


Z 


Y 


Z 


OM1 


7.91el0 






0.774 


0.008 


0.2787 


0.0220 


0.2503 


0.0207 


OM2 


7.65el0 


0.036 


0.0033 


0.782 


0.008 


0.2878 


0.0260 


0.2802 


0.0257 


AC1 


8.08el0 






0.767 


0.008 


0.2698 


0.0226 


0.2440 


0.0210 


AC2 


7.50el0 






0.780 


0.007 


0.2670 


0.0263 


0.2735 


0.0255 


OM3(a = 2.2) 


7.38el0 


0.0463 


0.0060 


0.755 


0.013 


0.2970 


0.0300 


0.2909 


0.0296 


OM4 


7.71el0 


0.031 


0.0022 


0.788 


0.006 


0.2820 


0.0235 


0.2724 


0.0230 
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Fig. 3. Echelle diagrams for the six models presented in Tables 2 and 3. They display the frequencies of the oscillation modes 
(ordinates) versus the same frequencies represented modulo the large separations (abscissae) 
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Fig. 4. Small separations for the six models presented in Tables 2 and 3 



0.11 and 0.19, and a mixing length parameter adjusted on so- 
lar models. Overmetallic models with a larger metallicity were 
inconsistent with the measured parameters (log g, log L/L Q ) ex- 
cept if we increased the mixing length parameter, which we did 
although it is highly unprobable: this is model OM3. Two mod- 
els (AC1 and AC2) were computed with the accretion hypoth- 
esis and external metallicities 0.11 and 0.19. Finally a special 
model (0M4) was computed with the precise external parame- 
ters of the Hyades stars. 

Observations of this star are planned with the HARPS spec- 
trometer in La Silla. Considering the precision of the results 
which were obtained for the star jj Arae (Bazot et al 2005 we 
are confident that it will be possible to distinguish between the 
different possible models from the asteroseismic observations. 



We will then obtain the mass, age, outer and hopefully internal 
metallicity of this star. 

These results will show whether t Hor has been formed with 
the Hyades or not. If it does not belong to the Hyades, the con- 
clusion will be interesting in itself. If the results show that the 
star belongs to the Hyades, it will be taken as an evidence that 
the overmetallicity is primordial and reflects the metallicity of 
the formation site of the star in the Galaxy. This will give a 
way to identify the right scenario for the exoplanet-host stars 
overmetallicity. 
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